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ABSTRACT 





Cross sections and kinematic distributions for the trident 
production process J= + z - ^ + c ‘ t e + t z {with Z = p and Fe) 

are given for beam energies of 100-100 GeV at fixed (e"e + ) masses 
from S to IS GeV. This process is interesting as a test of quantum 
electrodynamics at high energies, and in particular as a test of the 
form of the photon propagator at large timelike (four- momentum) 2 , 

B or this purpose, it is desirable to impose kinematic cuts that favor 
those Bethe- He i tier graphs which contain a timelike photon propa- 
gator. It is found that there are substantia! differences between the 
kinematic distributions for the full Bethe- Heitler matrix clement 
(which involves four graphs) and the distributions for the two time - 
like- photon graphs alone; these differences can be exploited in the- 
selection of appropriate kinematic cuts. The competing virtual 
Compton process (for Z - p) has been estimated on the basis of a 
Simple model; its cross section is at least two orders of magnitude 
smaller than the Bethe- Heitler cress sections of interest, at the 
energies and (e e ) masses considered. 




In .Ms paper cross sections and distributions for the trident pro- 
cess 

(j ^ + Z -* |j^ -t- e + e + 2 (LI) 

are calculated to lowest order of the conventional quantum electrodynamics 
for the esses of clastic scattering off a proton and coherent scattering off 
a spin-0 nucleus. This process, as well as the essentially identical 
reaction 

e* + Z - e* + p + p + ^ 2 (L 2) 

and the related process 

l* + 2 ~ +■ t" + i + + z (1. 3) 

(where £ = p ore) are of interest because they are sensitive to possible 
modifications in lepton electrodynamics at high energies. In particular, 
we ■ .interested in facilitating the detection of possible deviation from 
conventional QED. especially in connection with the mass spectrum of the 
lepton pair produced in processes (1,1) or {1.21. Such a deviation could 
be due, for example, to a modification of the photon propagator along the 
lines of the Lee- Wick ".heavy photon pole* 1 model. ^ 

In the Lee- Wick theory, the possibility of replacing the conven- 
tional photon field amplitude A by a complex amplitude A + iB u , 

^ JJ M 

where is anti- Hermltlan, is studied. Such a modification (coupled 

‘•v-L - ^ ' \ n '. i J3 " fc ***r . + TT j v i. 



with the introduction of a massive indefinite- metric fermion field) yields 
finite results for calculations of observable quantities in hadron as well 
as lepton electrodynamics. Total cross sections for the resonant process 


,(,* + Z - £* + Z / * 


£ 


(U 4) 


leptons or hadrons 


where y or e ( and B is the massive indefinite^metric boson associated 


with the amplitude B » have been given in previous papers 


3*4 


Hh 



The Bethe-Heitier diagrams for process (1. 1), in which only one 
photon is exchanged with the nucleon or nucleus, are given in Fig. 1 
(for Incident p~), There are four such graphs, corresponding to the in* 
teraction with the nucleon or nucleus Z occurring at any of the positions 
(a) -(d). The competing virtual Compton graph is indicated generally in 
Fig* Z(a), The virtual Compton cross section for scattering off protons 
is calculated in this paper using a simple model of the hadronic inter* 
action, It is found to be two to four orders of magnitude less than the 
corresponding Bethe -Heitler cross sections of interest. 

The main effect of the Lee- Wick (or similar) modification upon 

processes (1. I) and (L 2) Is to enhance the singly-differential cross 

section ^ in the vicinity of the resonance, where v is the (mass)^ 
dv 

of the timelike virtual photon. The enhancement of that part of the cross 
section which involves only diagrams (a) and (b) of Fig. I (L e, , the 
graphs which contain a timelike photon propagator), is given approxi- 
mately by the factor* 1 2 m (tm ^ provided v is not ^ rn 

o B r* 




1. 
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rMi ’•K. - C . r “Hi . >>• • ,•'*,»• tr« ■ ■»• » . " 

, - ‘ C * 

Thus if pft - ZO GeV, for example, &ft 80 % enhancement should be 
seen It V “v" - Ifi GeV. 

However, if erne mcasutes “ integrated overall kinematic con- 
figu rations consistent with fixed vj, it will be Been that this enhancement 

f£ . • " *■ 

effect will generally be nearly totally obscured by the much larger 
contribution from those diagrams [see Fig, I, cases fc) and fd)] in which 
all photon propagators are space like. The reason for this is simply 
that if no additional kinematic constraints are imposed, the sp.iceljke 
* piiotc.n joining leptonlc vertices in cases (£) and (d> oF Fig, 1 can be 

much closer to the light cone lhan can the timelike photon in cases (a,) 
and (b), for J v in the region of interest 10 GeVl* Consequently, 

V?e Bhail^aUo be concerned with the kinematic distributions tor pro- 
cesses (If 1) and (1. 2) * and with the types of kinematic cuLs that can 
be made to increase the relative contribution of the timeiike graphs. 

This information is useful as a guide for experimental searches for 
possible deviations of the Urn el ike photon propagator from the pre- 
dictions of conventional QED, The distribution^ are obtained effi- 
ciently by a binning procedure incorporated within a modified Monte 
*f Carlo integration scheme* 

In Section II, the matrix element (corresponding to the diagrams 

of Fig* t) is given for process (l. 1 ) [or (1.2,1 with the exchange ,4-— e]* 



Section HI gives the differential Cross section expressed in a convenient 
set of variable-. Section IV describes the phase* space integration and 
numerical methods used: In Section V the contribution of the virtual 




Compton graphs is discussed. In Suction Vt the rcstilts are presented and 
discussed. In Appendix A r expression^ for the relevant dot product it utsd 
related quantities are derived in terms 01 the variables defined in Section 
I1T Appendix 0 gives the limit of integration for each variable. ,nd de- 
scribes the method ared for incorporating kiflenruUc cum tfeo the integra- 
tion achemcl 


v, " / 
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U. MATRIX ELEMENT 


A, Lqptbruc Part 


Tiu! matrix demerit corresponding to the four BethD-Heitlcr dia- 
grams of JTig, i Ik 5 


U J = K V 
0 c 


K * r + S ; 

? s a 


(2.1) 

(2.2> 


where 


3-1 

T = ie v * u. "Y v 
jj „ 2 JU 3 


■ :, t a + Zk J 


* u b (V 4 *’ + Zk t o )’ r )J U 


(2. 3) 


is the contribution of the timelike" diagrams (a) and fb) of Fig, 1, and 




-ie^D * tt T u If [d “*(y q.Y4 2k. K 
l T 2l c W 2 a / T X 


V Zk 3 0 )] V 3 


(2.4) 


is the contribution of the H e pace like" diagrams (c) and fdj* Here k, k^, 
k^ arc the respective four -momenta of the incoming and outgoing ^ , 



and the e and e (see Fig, 11 T satisfying k^=k^=-m ^ and k ^ - 

1 *r 2 

k-, ~ ^ ^;Euid q - p 1 - p where p and p' are the respective initial and 

2 2 2 

final four- momenta of the nucleon or nucleus Z, with p - p' = -M . 


t 

.» , 


The propagator denominators are 


= -<W 


D = (k - k { )‘ 


- q - 2k*q 


= q + 2k^q 


D c - q + ?k 2 -q 


D d q + 2k^q . 


(2. 5} 


The spinors corresponding to the particles of momentum k, and k^ ( 

and to the anti particle of momentum are denoted u, Uj, and 

respectively; they are no rma I ^ed ao that uu - 2m and - 

- v„v_ -2m. 

3 3 g 


B. Hadronic Part 


The hadronic vertex and photon propagator q in Fig, 1 contribute the 


factor V „ The V used here are identical to the corresponding factors of 
C O 

6 

paper I. For elastic scattering off a free proton, the dipole fit is suffi- 
ciently accurate for our purposes: 


!, 





* 'if* ' 




;• >r~ 

Slnl 


whe 


and 


Then 


with 


■Bfij 


Str- ' 



V f F lY(j 4 \iy l ( Vb . Va) ; e , Up , (2 . 8) 


where % “ nd V * rc ,he -«* proton (spinors , satiB r ying 

Vp s V u p' = £m - 

For coherent scattering off a spin-0 nucleus of charge Z, we have 

V = (eZ/q 2 ) (pi p ') r ( q 2 j . 

0 0 

We Minnie for the nuclru. a termi charge distribution in the Breit frame 7 

4 (% - ■» 


P r W“ [ ! + exp [{r - r Q )/c] | * J 


r 0 *(l.)8A - 0. 48> ferml. c = 0. S5 fe rmi, 


12.11) 




and Cp normalized such that i ^p(r)d"^r = l. The form factor is given 
by (see paper 1) 


F(q 2 l = [ j a F (r)e‘ q ' r d 3 r) 6 (O.ZS GeV 2 - q 2 ) 


(Z. 1Z1 


integrated in the Breit frame, where 9(x) - 1 for x a O and 0(x) = 0 for 

8 

x <0* We have thus imposed the arbitrary but convenient cutoff 
q £ 0, 25 GeV Z . 


5 


For each of the above choices for V , we define 

a 


H = S 7E h ' V V ® 

a t * b 0 - 


u* n> 


where £ denotes the average over initial and summation over final 

spins (if any), [For any four-vector a = (a t ia Q ), we define 

-»* v . - • 

(a , U It where $ denotes c ample* conjugation, 1 Since H = H 

CJT T0 

and a H - 0, we may write 
ff at 


„ I./, j siA.hL 

" 1 ” ,*/ “'V 


-=-/ D --^5- 


« z %ll P T Z ■% 


*: E^a. 


Thus, for elastic scattering off protons, we have 


H 1 = ^( F i + 2 F z ) 2 

H z = 4 ( f i m2 * F *V) : 


12. l?l 


V I 

j. (2.141 
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and for coherent scattering off iron. 



H z = (ZMZFr 







(Z, 16) 




111, DIFFERENTIAL CROSS SECTION 


The differential cross section for process ( l. 11, assuming the 

beam to be unpolarized, and summing over spins of final- state parti- 

5 

cles, is given by 

8 d 3 p * d 3 ^ d 3 k 2 d 3 kj 

4? ^l^UbVVWSo 

M * f* 

% *, ^ * ’ V*-. rr-rv r 

' x s (4, tit - k 1 ^<y- kj -'qj r/iV '<sif ? 


where denotes the average over initial and summation over final 

lepton spins. To simplify the six-fold numerical Integrations fd:& 

& 

cussed in Section IV), it will prove convenient to express d .r in 
the following set of variables. Wc define 


s k*q 

= 2tk - k^l* q 

u } - (k - kj) * tk , - kj> 



i" ; ' ’ - .T|C. * i -5 

0 / = azimuthal angle from an arbitrary axis io p ' about polar axis k 

P 

in the lab frame (p = G), d ^ ^ = azimuthal angle from p to ^ about polar 
axis k in the frame k = q (called frame “A"), and $ azimuthal angle 

from k to k^ about polar axis q in the frame £ - q (called frame °B"). 




i 

■ 

«• 4^ 


> M Ct. 


dt <!u l *t5' 0 7 

pvd U ,d S , A 


\l*$l 




V TtnS^ %£ 

* 1 y* 


* 

L? - , "TG W- ,\. L ^ '*- 


TJ ~7_ ** v 1 



du _ d# jM 


'l !1 ' 




-rV 


,/* 


t 2 




- 




where 


4, A 


*W 


i ■ - • 


- i © 

t / K K W , 




•t c.4Sf! ; r^jj fSa&ifc*; 

; r v J? - 

l5*>k r 


' *0. *T 


£l«> 

it 'V tV '> „ ^ 

(3. 5) 


^ ' u >> ijfc $fr.- * 

#rf«t * , :*’-•, 


i . . 


c ,*^t 




and ^cri: thi* iwr braetust* are SfSt'ft respectively, to 


• j^V P M Jt ]■ 

k p p 0 -ta>, ’ :t -l^lOM ' 









’ * a nicl 


i ;; ' V VS-'^VSou ■ 


'M'Mp'fa W-~S ~ . Sifii . i 


' ' - ' I B| 2 ] . 


[i4ntt}‘ 4 -jt - Lb 

#^®HW|F "W*:- 


>aia • 


q 1 * - q L - 0, 15 q, (2. 2 4} implies 

3 , ax T ct 




"2 

H sH.li f— | p L p 

aT o 1 1 oq ^2 (, Qt- t 


f 3._6v 


j4^y 
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;- 3} ■’ 

■ J3S«?Z~ , y.j9B , ^r Vnr 


L ^ and p L, p are evaluated in invariant form, using Veltman's 

9 

symbolic maniputation program SG.HOONSCHIF to perform the 


^ tedious truer calculations, (It should be noted that the expression 

for L is not needed m the ease of coherent scattering, since H. = 0. ) 

rr ■_ 


The resulting expressions foi 1 L and p L are composed of dot 


products Involving the |our-vcctdrs k, k^ f k^ t q, ^nd p (the latter 
appear* only in p L T P^)- These dot products are given in Appendix 


A as Rmcticms of tfu ^eycn variables v,t, u u^ % ^ and In 


the next section v *■ U mad e **? the fact that only kj* p, fc^'p> and k^-p 
depend upon d ] A - *^nd only k^-p, k^-p, k«k,, %'S k^, k-k^, and k- 
depend upoii p ^ The confutation time is enormously reduced by 

t Md Hz 1 


faetorinE; out th« and $ dependence. 




' 




Jikic a? 














• -4- 




' * u ^"- 


^j0L 


I 



!a Ibis section the procedure used to integrate the differential 
Gt&ss section over -alt variables, (except v, the square of the invariant 
mass of the electron pair] is described, and the lab variables over 
which binning is performed are given. 


A* Choice of Integ rati on Variables 


The kinematic hr: jits for the variables defined thus far (L e. , 

L t Uj* ti^, Uy and arc given in Appendix B, Two primary 

advantages have been gained by introducing these variables, rather 
than attempting to integrate the differential cross section in the lab 
frame, f irst, the fact that th^ integrand depends in a very simple 
way upon and A enables a considerable reduction in computing 
time, as mentioned in the previous section. Second, if the integra- 
tions were performed in the lab frame, special procedures would be 
required for choosing the points of evaluation in such a way as to take 
account of the relativistic peaking of the integrand in the forward 
direction, This difficulty is eliminated by doing several of the integra- 
tions in the CM frames of selected subsets of particles. (Frame M A ,r 
is the CM frame of k ,, k , and k ; frame ,r B", of k and 1 , ) 

j M i 4 J 



However, the integrand is still not sufficiently smooth in these 
variables to perform a Monte Carlo integration of sufficient accuracy. 
This is because several of the propagator denominators can attain very 


15 





small values in certain kinematic regions. It should he noted that each of 
the denominator a [D, D , D , D , D , and l) is linear in v t t, u , 11 , and 
Uy and is independent of and 0^ [see Eq. {A, 15)]- This fact helps to 
make it convenient to transform to a new set of integration variables, chosen 
so that the new integrand (including the Jacobian) is moderately smooth in 
the new variables. This is of course equivalent to optimizing the placement 
of the points of evaluation of the integrand, to take account of rapid variations 
in the original integrand (resulting both from the propagator denominators 
and the hadronic form factors). 

The following transformation has proved useful in this regard, although 
a somewhat different choice could also have been made: 


P G(tO dt', for scattering off protons 
t 

1/7 

[I* (0) - I y (t)] , for scattering off nucleus. 



14 . 1 ) 
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G is the dipole form factor [£q. (2, 7)], and F is the Fermi form factor 
[Eq- (2, 121]* In the case of coherent scattering, V is tabulated as a function 
of i, and the conversion between t and l 1 is accomplished by interpolating in 
this (sufficiently dense) table. (See paper I and Ref* 7.) 

The rationale of the transformation (4* 1) is as follows. The Jacobian 

fur the mapping l-*Xj reduces the peaking at small t due to the hadronic 

interaction factor V * Next* 

1 


d* 3 = (2uj - v)(DD h Y l dii 2 


(4, l) 


results in the complete elimination of the D~ 1 and D, * factors in the cross 

b 

terms between diagram (b) and either of diagrams ic) and (d)* in Fig. I. 

The effect of this transformation upon the entire factor b H is more 

0T ctt 

complicated, but a reasonable smoothing of the integrand results* The 
**■ je^ transformation, whose Jacobian i* given by 


d %= D c D /6 


£*• du^ 


(4* 3) 


similarly smooths fin partj the u_ variation due to D and D .* The denomi- 

> c d 

natora v and require no special treatment* since | | > v and we are 
interested in the case of large v. 


B- Numerical Method 

We obtain the singly- differential cross section do/dv* and its distri- 
butions in lab energy* angle and related variables, by a modified Monte 
Carlo method. The six- dimens lonai Integration region ( in the 


17 


variables x, , and ft , with the outermost integration being 

l-*4 3B 1 A 

x^ and the innermost, £^1 l® partitioned into 

N = N. 
i=l 1 


14 . 4 ) 


rectangular boxes, corresponding to bb equalLy-spaced intervals 

along the x. axis (j*I-4), N $ along <£ 3B * and along ^ For each box 

in x ^ * 2 x^ x^-space. two points (x^, x^* Xy aral f x |* *3" *4^11 

are chosen pseudo- randomly* The bulk of the evaluation of and 

p L p . which consists of several thousand operations, is carried <*it at 
r fj ot t 

this level of the integration. This is possible because the relatively 
simple dependence upon £3^ and $j A can factored out* 

At each of the points chosen above all expressions independent of 
0 ^ and are computed, N^ pseudo- random values of $ ^ ar ® 

selected (one in each of die interval®). For each value* k- ky k \ k 2* 
P- k , and k^k^ are calculated (see Appendix A), and all expressions 
independent of are computed. Finally, pseudo- random values of 
^IA ^ 0Iie 1>ef ' nterva ^ are selected, the integrand is computed at each 


"1A 


valuf , and the values of q and (In the lab frame) of kinetic energy 


Ty log 10 (1 - cob 0.) where 0^ is the angle made with the forward direc- 
tion, the transverse momentum k ^ = j k | sin 6.* and the quantity A k > * 
(k| (t - cos g .)* each for ! * 1, 2, 3, are computed* (See Appendix A for 
details of the T. and cos 0, computations*) The differential cross section 
for each point (I* e* , each simulated "event 11 ) is assigned to the appropriate 
one of a large number (>. 10000) of internally-carried bins in each of 


IF 


these variables* At the end of the calculations histograms are generated 

in these variables. Each histogram contains about 25 bins, which are 

- 4 VI Hfeo 

selected so that an approstimately P <j ua l partial cross section :iQp) 

falls within each bin. 


Because the vast majority of the inkgi and caiculatioj^ are done 
inside only four of the six integrations, the bulk of the computing time 
is spent in the calculation of the bin assignment* rather than of the 
mtegrano. For given values of v and inoummg muon energy £, the cal- 
culation of — and Us? distribution ir. and m the iant .atory variables 

given above* requires about a minute of computing on u» l&M JGa/95. 

Lpr • ^ *■ 

MMH - ' yrlF . ^ . >:r ■ ; ■ *v 

Since each set of * J-4 values is selected lindfcp* ndviitly an,} 

randomly, a statistically unbiased esLimata ofthe 'accurac y "at the :< • 

• . VT- ' . ‘ , l \ 4 

integrations is given by the variance 

(o.d .) 2 . *i. 

*S$ ^ ,>££'. w ; ; \* ¥ . .. 

summed over *11 points in '«!»«} where I Is the Integrand of 

the x l-4 inte S rat “^8 (already integrated over © iQ and and "s. d. - 

denotes the standard deviation. The total integral is. given by the 
approximation 

~ %, uNjtyijjy- 1 rf(^,„ 4 > 3 ♦ . 14.6} 

For the results presented. <s.d, >/( Mtfv\ is approximately 2% or less. 



The reproducibility of the values of dg/dv and of the kinematic 
distributions He. , for different sets qf pseudo- V and pm numbers) has ‘ 
been tested for -tcv^raf cases; Ihe dc/clv values given are reproducible 
to within about 2%. and each of the approximately 25 lun population* ffor 
each distribution) to within 10%. An independent estimate of the accuracy 

of some of the kinem^uc difitri butioa> lx obtained by noting the extent of 

■+ * ' 

the discrepancies between the v and > distributions. which by symmetry 
would be identical in an exact c dlcblation. The differences in the e* and 
e~ bin populations are found to be on the order of t several percent. 


y.1 1 1,' 


,;k 








te- w*- 


i *rj2 ^ :j 








' ; t jjU'i 

^VW7V ip77 >T '~ 1 V vc 
: 1 V^J jWgC^M« ■;>' -^Vf/' 

~*-V.. * 
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V. VIRTUAL COMPTON GRAPHS 


th-*n dotfs Uc rtag.in detail the complexities »f the Hadronic inter- 




cuon in ihoAwf>-phou»r>- rxcharirge tftaft [fig- 4Wr 1 wfc have calculated the 

out ri but ion of the virtu ai Compton graphs for scattering off protons* according 

o a simple mode t . * J In this model. tfi> proton is treated as though it were 

are, except that an overall dipole form factor G[(p -pi*'] [see Eq. (2-7)] is 

2| «r ^ajpao A# % ■ ^ u. - -Lcf**”* r"“L ',i . ■- 4**“ ^ r ^ • 1 '• v 

.. icorpqrated into the matrix element. Apart from the factor G, the tws virtual 

ompthli^aphs are then given in Figs* 2(b) and (cf, Under the substitutions ^ 


m 

u 

k 


fcftg* 


* 

P 

P' 


(&- 1> 


iesc diagram ; become identical to graphs iaj-and lb) ol Fig, I (for Z - bare 
r ut on) . Th!i ci roumstancc makes It pass iblt To calculate the singly-dlfferent ial 

rose section for virtual Compton scattering by (he same method as that used in 
rctlbns II and ill (and the Appendices}, if one replaces by in Eq, (2. 1). 
■ekes the substitutions (5. \). and includes the appropriate form factor. The 
*ri» format ion of Section IV is replaced by a more appropriate one 


Xj = *it 


£ 1 

>,■/; 


X 4 3 U l 


G (o^ # - v - ij du^ 1 


15.2) 

lit la 




IHHB 1 




KAf- •'>*, ■- - », vs^T*-- 



where t and u i (i * 1, 2, 3) are defined as in Eq. f'3. 1) with the substitutions 
(5, 1) applied [thus t = {kj-10 2 . etc.J, The v.rtual Compton graphs will be 
found (sve Section VII to make a negligible contribution to the cross section for 
the large v values of interest. 


u 








VI, RESULTS AND DISCUSSION 


Wv have calculated (do/dv) and its distributions in q and the Lab 

quantities T,, cos iL, k and jk_ f for several types of cross sections 

g of interest, « and - denote the cross sections for the Bcthe-Heitler 

p Ft 

process [Fig. 1(a) -{d}] for clastic scattering off a free proton and coherent 
scattering off iron, respectively. The cross sections corresponding to 
the two Beth e-Bei tier graphs which contain a tlmelike photon of mass Jv 
[i. e. , graphs (a) and {b) of F ig. 1, which are calculated by replacing K 

0 . 

by T in Eq, (2, I)] are denoted by g * and a _ 1 . The Cross section for 
o P re 

virtual Compton scattering off a free proton (calculated according to 
graphs (b) and (cl of Fig, 2, with the dipole form factor G[(p # - p) 2 ] 
included in the matrix element) is denoted by g 

P 

In Fig. 3, dWdv and d.-^Vdv are plotted versus v, the (e% + ) 
mass sq larad/ for incident muon energies of 200 arid 300 GeV. 


As discussed in the fntTodv ctlon, we are interesteci in the use of 

/do ,/da \ 

kinematic cuts to increase the ratio ( 1 /( ~g^~) ( an ^ similarly for 

Z = Fc), It is efficient for these purposes to make such a cut in coa 0 ^ 

where 0^ ^ is the lab angle of the outgoing muon relative to the incoming 

beam. For reasons to be discussed below, the somewhat arbitrary cut 

cos 0, , . s 0.%6 has been chosen for further study. We define a a a 

1, lab p ( c 

that portion of j which arises from aU kinematic configurations satisfy- 
ing cos o , , i 0, 996; o ' i t- i and au are defined mutatie 
9 Ulab p P c Fe*c ^Fe , c 


23 


For electron pair masses of /v - 5* 10, and 15 GeV and incoming 
muon energy E ^ 100* 200* and 300 GeV, Table I gives tl/ZHdg/dv) 

where g denotes the five types of cross sections cr^. cj^i and 

2 

y* . The average values of q and the Lab quantities cos 9^, 

cos 0, k . . k * r 4k ni . and ^k.. are also given. { The subscript 
2, 3 1 1 i£, y |j 1 3|2. 3 

lr 2, 3” indicates that the electron and positron distributions have been 

averaged together; this is reasonable since these distributions wcuid be 

identical in an exact calculation. ) For three of the (/v, E) pairs. Table I 

also gives ( 1/Z)(d g/dv) and the average values of the kinematic quantities 

for the case of a , n * . , and g * 

p,c p, e Fe, c Fe,c 

For the case /v = 10 GeV, E = 300 GeV, Fig. 4 gives {in histogram 

form) the computed distributions in and ^ for each of the singly- 

differential cross sections dg /dv. dg f /dv, dg /dv* and dg ' /dv. 

p p p.c p,c 

Figure 5 gives the corresponding distributions in log^U-eos 8^) and 
log lo (l.c°s0 2 3 L Fig. 6, in and ^ and Fig, 7, in q Z . 

We list several results of the calculations below. 

1, For the {/v, E) values of Tabte I, dg /dv is at least an order of 

P 

magnitude greater than { l/26)ida /dv) [ and likewise dg 7dv >> {1/26) 

ft p 

(dg * /dv)] except far from threshold (at /v = 5 GeV, E> 200 Gev). 
r e 

2. An essential check on the calculation of dg'/dv for 2 - p and 
Fe is obtained by comparing do'/dv with the total cross section g(B^) 
for production of bosons of mass /v via process (1. 4L By comparing 
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the results of Sections 11 and 111 with the cross section formulae Of 
paper I t one finds 

a(B°) = (1-rvfa )( 1 + 2m Vv|’ 1 (I - 4m^/v f (doVdv) * 


(5. 1) 


The numerical results for d?'/dv satisfy Eq* [$* 1) to within approxi- 
mately 2%. 

3* In gene rail da /dvi IQQfdo' /dv) > /dv). This can be 

P P P 

roughly understood in terms of the propagator denominators. In 

d^p/dv, the "spacelikc" term S [Eq. (2. 4)], with Us overall factor of 

D ~ (k - in the denominator, dominates the matrix element; while 

the matrix element for da' /dv contains the denominator v instead. 

P 

Ignoring muon mass terms, D & ZEAk^, which typically is several 
GeV* for da^/dv (see Table J), Therefore one would crudely expect 

that (da /dv}/(da' f dv) 2:0(v 2 /D 2 ) s= Of IQ 3 ); in reality the ratio of 

P P 

cross sections is closer to 100. 


To see why the virtual Compton cross section da* /dv is so small, 

P 

we consider the factors arising from the photon propagators and the 
hadronic form factor G(q Z ). For doWdv these factors give v V^G^); 
for da^/ dv they give v ^G^fq 2 ). Using the average values of these 

quantities (obtained from Table I)* one finds (da' /dv)/(dg* / dv) == 0(100). 

P P 

4. For the case of do^/dv and dj^/dvi the outgoing muon typically 
comes out nearly forward and has a moderately flat energy spectrum; 
white the electron and positron come out at large angles from the forward 



mmm 

■MMH 
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direction and have spectra which ‘are strongly peaked at small energies: 
(See Figs* -1 and ?* ) The value p£ (T, ^ ;s somewhat misleading 

in this respect, because of the highly skewed distribution, •Mtho-sgh (for" 
do^/dv with /v - IQ GeV, E - 3G0 GeVj the mean 3 is } 9 CeV and the 
median is ^ 85 GeV, 40% of the ^distribution falls bek-w 15 GeV. A 
similar situation occurs for other (/v ♦ E) valuer 

5. By way of contrast, dr '/dv (and d-j' fdvl exhibit distributions 
p re 

very different from tfTdle above. "fhe outgoing smmn i 4 ? or. th,* aver 
slow and comes out at- a large 2U f; 

are fast, have moderately flat spectra and corf ' <*< ut at a siWm whit 


Icctron and |fc<Hi*t>n> * 


■ 

: 4^. * 'V iV-'T, v ’Qp? i . "* ' '>• TV" ’ 

angle (er 5* ). This is the situation wluch would ^rcur in the c^se cf ro$tmantj3 
electron pair production, via process (1.4). *1 


6. Because there Ib such S slight overlap between the iog 


10 


- c os « 


distributions of da /dv and da '/dv [see F ig. a>J. one can choose •* cutoff 
P P 

. 

value of cos 9 . , called c '. such that >98% of d? /dv ’its ibwTc'; 

p - .? 

>50% of dcr '/dv lies below cV A suitable (though not uttiqufe) choicr. * 

P 

based upon Fig. 5{a) and similar distributions at other (/V E)" v.uucs 

is c* c 0.996 (for which Logoff - st^.4). For this cutoff. Inspection’ 

of Figs. 4 - 7 indicates that the da /dv distributions are generally 

p* c 

similar to those of da Vdv ■ while retaining seme of the do /dv peaks that 
P p l? 

arise from propagator denccnlnat ur* . 

^ r4(J ■ ’* v * •- v V' v * jaMBr-CJ* 




*> 

C55 


because of the large ijuaUlative difference* between the drr /dv and 

P 

4 ‘ 7 p c. nchj^e that It is possible, by Impbsing the 

*pp r « r rUHv ktttW>»afie cut- • tifcg*»ted by F.ga. 4 - ?, to.deslgn experiments 




v/bich are highly sensitive i, possible modifications in the timeiike-phaton 
Bathe- Heitlec tfift* aViar . e value* V v, ’ the (o' £ ) mass squared. The 
Competing contr i button ^h«-' epacclikc^hotori graphs can be made com- 

parable to lh<- tifrielike contribution. «whU e the virtual fjspmto* cro*» 


section l* extremely -sipalj. according to the simple model used here. 

■ . - 

Althiv^h a cut in cos v has been considered f«y- purpose* of Illustration. 
*# 




«V- 
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APPENDIX A 


Using the definitions of Eqs. (2. 5), (3. 2) f£, , and {3. 5}, the dot products 
of k, k^, q t and p are calculated as follows. We have 


k-kj = 


k*q = 

u i 

k*p = 

- Mk 0.1a b 

V* = 

u i'z u z 

v k i 

Z 1 

- m ~ -r \ 

e Z 

N_ 

II 

1 1 i 
'I t+ 4 U 


.)- 


. i ■ i i 

k 3 .q = - 2 t+ 4^^“j 


q-P = ' S'* 


(A. 1) 


In frame "B" the polar and azimuthal coordinates are denoted 9 and <i \ 

o Q 

the z axis is defined to tie along the 9 n ~ 0 ray, and the x and y axes satisfy 

a 

(@ B * and |~jj . ^), respectively. By definition, k yB = <l xB = 

3 0| WC the sense of such that a 0. Thus 


k* k* = k k_ + k „k - k_ _ k_ _ „ , 

3 xB 3xB ZB 3eB OB 30B 


(A, Z\ 


We have 


where 


thus 


and 

can be calculated. 


and 


one obtains 


k 


q GB = & ( fc ' i u z) 




k„ n = JL 

0B 2/7 


T| = v 4 t - u 2 + 2n % ; 


k zB * ( u I*\b%b)' >Mb 


k *B = ( k 0B ' k *B-* 


■’) 

u * 


m 


Using 


k 303 = 2^ 


k 3zB = ~ V/ft ' 


Ml 




3B 




(A. 2) 
(A, 4J 


(A. 5) 


(A. 6) 


fA, 7) 






" 

, , , - * - j| 
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Therefore 


k-k 3 s 


4 + 


\V'+-<t} (; 44 ) 


1/2 


cos «. 


3B 


where 


C = *«,♦ n(*4 u a ) • 


(A. 8) 


(A. 9) 


The other dot products that depend upon 0 but not upon p 1 A , are 


"1A' 


k ' k 2 = .-r, - k-k 3 
V '2 ! U 1 + 4 u 2 * 2 w 3 ‘ k k 3 


and 


v k s = i( v + * +u 3 ) -i u 


+ k' k. 
2 3 


(A, 10) 


In frame "A" the angular coordinates are denoted 9 and a ; the 

A -A' 

< # y, and z axes for frame "A" are defined in terms of S and * as 

A A 

fe scribed above (for frame By definition, p = k = k * 0 

yA xA yA * 

nd the sense of is arbitrary. We have 


P SxA P xA + k lzA P zA ’ ^ 10 A P 0A 


(A- ID 
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where 


with 


end 


Also 


k ioA ( u i + m u B -i v - i*) //8 
_ [(V*)^? 4 ^ ( t+ U 1~ U 2 )] 


IzA 


IxA 


2(9C) 


1/2 


(v V - H 


2 2 
- m v 
u 


IzA 


cos 6 


IA 


P 0A " ( Mk 0,iab ^ 2 l ) 1 ^ 

p *a 3 [ Mk o, i,b(‘ ■ u i ) ‘ \ *( u i + m u 2 ) / iB c) * 
(MW- Mk 0.1»b fa l ~ 4 lm / tZ ) 


M 


1/2 


(A, 12) 


6 = 2 u, 4 m - t 
1 y 


Z ^ 2 

c = u + m t 
1 u 


(A. 131 


V P = k 3xB p xB + k JyB P )rB + S*B P *B * k 30B P 0B 


I 







. A. '• tr • ‘t* ' * 

[> t ^‘iT,,^^’'^ y*J. 


where 




Poo s (^o.ub-r • S-p) '^ 

* - 



P in * [( il ' u 2 ) p ou • 1 " 






P .nB ( k Lb t> 0ii (..lab ' /.>' '/■ \)j *a 


wm 


I l .,2 2 2\l/2 ... 

PyB ( p 0B * M ‘ P <B * P*n ) ,A H ’ 


nd k_. 0 , k_ ire given by Eqa. {A. 3> through (A. ?}♦ i&r * - G. *,**'«.- 

3 A d Afl 


7^4Ss 


The propagator denominator c are therefore [nee E«i> f*. W] 

T>rf« ;vj«saR: «*rr»- •--> . - 


P ’« Ug * V - t 




D 5 t * 2u. 

*>-/♦ 2uj- 




D « ~ II 4 U 
c 2 2 3 




°d = 


IX. V*l 


Using these mulls, it is straightforward la derive express ions for the 
b variables oyer which the kinematic distributions are obtained, The 
aetlc energies of the outgoing Leptons are 




m H9>*rC 











$S§ M' 

“ Jr* V -i 'T' ^ w 'J^ 1 * ,r 

— •' *33 



« 


'"V "V ’ 

^..V. ;ll 'M <i=l. 4.31 


Also 



"M 


cor-J 


k*k t 4 EE. 
i I 


T - 


r k i; v, 


, , - 


. '^Vi 


where E i *^*£>4 . 1* ifc* '% f " 




,<W*. ~ -XL 




fA. m 
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The integration durnains for the azimuthal angles V . C 3B , and 4 U . 

cachVCuaJ from ii :c 2” . T>« . , integration. is trivial. Since all computed 

P 

quantities are invariant und«*r the mapping 4^— x ^ ^ convenient to 

7 

integrate over ? ^ A Irurn 0 to *“ only, and multiply by two. 


The limits of th? integration arise from the condition 




|co. 9 jb | ' l 


ib. n 


where % “ (q , Ij in frame "IT (in which frame k - y ql- One obtains 

3 Jf5 3 


** fyi <► 


U i- 6 *3 £ 11 34 


where [see Eq* (3. 5l) 


V- 


1/2 


IB. 21 


Similarly, one require e 


I c °» « 1A I * 1 


o. 3i 


where 5. s L'C* . k ,1 in frame "A 1 " (in which frame k = ql. Therefore 
IA l 


V <u 2 * V 24 


Vt *' ; 


: ir? . J: X ; 

ri*55 sa^-a yi 




' '* ••• ■ -. ... . c* 
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where [see Eqs, (A, 13)] 

I 


u 2±' t+ U 1 * j( v ' rn J Z )( u l- 1 ) 4 (atij-v-u) 2 . 4mj 

The domain of u ^ is given by 


m i / 

a/. 


IB. 4) 


“i- S “i S U l4 


where Uj _ is the value at which the physiol u, region vanishes. 


. -ivti 

1- 2 V * 2 


IB. 5) 


and where u J+ corresponds to the nucleon or nucleus coming out forward in the 

lab frame, 

.11/2 


1 4 ft/zM r 


- k 0.1,b t/2M 


u l + - I* u [' 

Requiring £u J+ yields the kinematic limits on t : 


(B,b) 


t %t % t 


whore are the roots of 


at 4 bt + c = 0 


IB,?) 


vith 


1 

w + 

^0, Jab 

2 

m 

4 -JL 

4 

2M 

2 



4M 

/V 

, lab 

Vi \ 

l 

M + V 

\2 v4rn u /v ) 


(tv 4 m u a) 2 . 






f B* 8) 


* 



The limits on 11 ^. u^, and t are similar to the corresponding limits in 

production (Paper ll , since the mass (/v 1 of the virtual time like photon 
(in diagrams [a) and tb) of Fiji- 1] fixed. This renders the kinematics similar 
to that for the case of three outgoing particles (nucleon or nuc lens . lepton, and virtual 
photon), if one is not specifically interested in the k^ and distributions. 

The threshold beam energy for given v is 



i(i v * n \ /v ) 


4 m + fv 


l B.9) 


It is sometimes useful (see Section VI) to be able to integrate over only 
that portion cf kinematically available phase space which satisfies an arbitrary 
constraint such as 


cos 0. t < c 1 * 
1, lab 


(B. 10) 


An efficient means of doing this* which avoids calculating the Integrand where 
It Is unnecessary, is given below for the particular constraint (B. 10). Cuts in 
other variables can be treated similarly. 


Equations (A. II) —(A. 13) and (A. 16) - (A. 17) imply that 


cos 0 , , . is of the form 
li lab 


cos H 


V^_ t0 * # !A 


*>* e < e °**lA + C S”* *1A 


J7T 


1. lab 


<B. 11) 


l 
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m 


where the c, (I = I. 


5) are independent of i* lA and Therefore 


a ( co * -ia) 


* 0 


(B. 12 ) 


U and only if 


cos $ 


C 1 C 4* Zc E C 3 
W ' C 2 C < - Ze l e S 


IB. 13) 


For each choice of t r u^* u^, and u^. the values Of cos ^ obtained by setting 

C0 * ^IA to + 1, - 1, and the BHS of Eq. (B. 13) [if the absolute value of 

(B. 13) is < 1] are each calculated using (B. 111, if each of the resulting values of 

cos 8 j exceed* the cutoff c # of (B. i0) f then no possible * JA can satisfy 

IB. 10), and the entire phase space region corresponding to the given values of t 

»«* must fall outside the cut. If, oaths other hand, some # lA satisfies 

(B. 10), then the calculation of L and p L n proceeds as usual, and each 

00 o ot t 

pseudo- randomly selected is tested separately for the criterion (B. 10). 
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10 

This simple model was suggested by T, D, Lee m the basis or 
scaling arguments; the model is intended to be used only at the lar| 
masses considered here. 

11 

The formal substitutions (5. 1) are used only to simplify she descrip- 
tion of the virtual Compton cross section rompv *ations~ In Section VI, all 
variables revert to their original meanings; thus, for .:•* 1 >n v 

VI denotes the hadronic momentum transfer fp* pf, ami nn\ the n.uon 
momentum transfer (k^ - k) [which q would have represented under the sub* 5 
stitutions (S, 1)] . 

■ - a 1. ; r f ^. 1 /* — tU * ** “ ^*1 

iVcr au ^ .v4l 


>*% ■* 



i#jc 1 2* ? r 

' -h-jj 

1 


TABLE CAPTION 




eSafe* 


Tahiti I. Singly^ differential cross section per proton {1/ Z)(do/dv) and 

aver*#* of kinematic quantities, for process (1, I) with 

_ 4. 

Z = p an 1 Z Fv. [The quantity 1 v denotes the (e e ) pair 
n i a « j squared. ]' Several types O; cross sections 0 arc co^ 

sidercd; 1 and 3*. arc calculated using the four Ucthe- 

p Fc 

Mcitlvr L»raph.v jI ig. J; *' and 0p r refer to the two 

ike" Uethc Hr ;Ief graphs only [Fig, l(*i and (hi]; 





- 


e p.e*k?#.,p' ;.c' awl «F 0 ,c 

err , cUun. vuth tjje ctitaos 996 Imposed, 

denote* the virtual Gomuton f r<*s*« section (Fig- 2), 








NaJSfei ‘ :S 







Fig* J. 


Fig, 2* 


Fig, 3. 


Fig. 4. 


Fig. 5. 


Fig. 6. 


Fig, 7. 



FIGURE CAPTIONS 

Bcthe-Heitler graphs for process (1. 1). Tht virtual photon 
carrying momentum q can be attached at any of positions 

(a) * (dj . In diagrams (a) and (b), the unlabc led virtual photon 
is timelike; in (c) and (d), it is spacelike. 

Virtual Compton graph for process (!» 1>* fa) General structure. 

(b) and (c) Graphs under the assumption that Z = bare proton. 

Distribution in fe e + ) mass squared, v, for process (I- l) with 
Z - p and E - 200 and 300 GeV, for the four Beth* -Heftier (BH) 
graphs (- J and for the two ,, tfmeHke M BH graphs alone P* - *). 

Normalised kinematic distributions in lab kinetic energy for' 
process (1. 1) with Z - p. /v = 10 GeV, and E » 300 GeV, for 
the four Bethe-Hcitler (BH) graphs!—) and for the two "timk- 

like" BH graphs alone ( )* [The quantity v denotes the 

(e e + ) pair mass squared. 3 (ii) distribution with no kinematic 
cut imposed; (b) Tj ( with the cut cos 9^ *0.996; (c) T ? no 
cut; (d) 7 , with same cut. 

*1 * 

Same as Fig. 4, but for distributions in the lab-frame variables 
log 10 (1 - cos B } ) and iog l0 fl - cos 0 J* 

Same as Fig. 4. but for distribution* in the lab tram* vers* 

momenta k . and k . . * 
i* 4.2,3 

Same as Fig. 4, but for distributions In q" k 
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TABLE l 


/v 

IGaVJ 

c 

(GaV) 

Type of 

Q 

1 do 

£ dv 

T 1 

ictv; 

T Z, 3 
<G«Vt 

cot 9 J 

co. e 2 3 

(CeV) 

k iZ. 3 
fCeVj 

4 k ll i 
(GeVJ 

4k | : ;z.3 

(G*V) 

2 

q 

(CeV 2 ) 

s 

100 

0 

p 

3,85 

39, 1 

30, 3 

0.9985 

0. 769 

0. 38 

0.89 

0, 0058 

0. 120 

0. 322 



°F e 

0. 0486 

IS. 0 

41.0. 

0. 5986 

0. 702 

0,26 

0.72 

0. 0046 

0.079 

0, 051 



0 ' 
p 

0*0578 

8 . 1 

45.8 

0, 9304 

0, 987 

0.64 
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0.072 

0, 246 



t 

°F* 

0,00180 

3.8 

48. 1 
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0.0203 

0. 066 

0. 048 



a * 

P 

0. QQ0270 

80,8 

19 4 

0.9995 

0.941 

0.58 

1,92 

0. 0083 

0. !95 

0, 899 

1 

zoo 

° P 

14, 1 

100,1 

49,8 

0.9994 

0.819 

0.41 

0.84 

0.0031 

0, 083 

0.225 



°Fe 

8,00 

43.9 

78, 0 

0. 9996 

0.761 

0. 30 

0.68 

0. 0025 

0. 042 

0.016 




0, 159 

17.0 

91.4 

0.9549 

0* 996 

0.76 

1,93 

0,0443 

0.038 

0. 160 




o no 

9, 1 

95.4 

0. 9904 

0. 996 

0. 32 

L 86 

o. oils; 

0.033 

0,014 



# 

O. 000422 

155,4 

22.0 

0,9999 i 

0.935 

0,61 

1, 92 

0, 0044 

0. 186 

0.869 




0. 112 

4.2 

97. 8 

0. 9327 

0. 971 

0.61 

la 64 

0, 0623 

0.033 

0. 161 



°F<r, c 

0.0741 

2a 0 

98. 9 

0. 9798 

0. 900 

0,27 

1.38 
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3. 3 
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0. 55 
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O. 033 
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1. 512 



r 

< 
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M 

°P 
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10 
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% 
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0, 0069 
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0.0611 
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0, 9645 

0- 989 

0,34 

3.69 

0.0259 
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% 

2.68 * 10' 6 
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0.9995 

0.971 

0,73 

3.74 

0,0105 

0, 226 
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°p. C 

0. 00180 

6, 1 

96, 8 

0. 9629 

0.898 

0.81 

2. 56 

0. 0667 

0. 13 1 

0,453 



a Fa, c 

3. 08* 10" 6 

3. 3 

98.3 

0. 9676 

0 .736 

0.43 

2.ol 

0.0337 

0, 128 

0. 122 



v 

0.000622 

4. 2 

97. 7 

0.9175 

0. 989 

0.70 

3.71 

0, 0820 

0. 129 

0. 464 



*Fa.c 

9.71 x 10' 7 

1*8 

99, 1 

0.9234 

0, 989 
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3.66 

0.0373 
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99 . 0 

0.9991 

0.774 

0 * 55 
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0 , 0054 

0 . 140 

0.407 



<¥. 

0 . OOI 2 5 

54 1 

122 . 9 

0 . 9*88 

0 . 7 Z 6 

0 . 42 

1,30 

0. 0049 

0 , 104 

0 , 083 



v 

0 * 00 20 > 5 

ie,* 

* 40 . 7 

0 . 9460 

O, 994 

0 . 96 

3.76 

0.0596 

0 . 092 

0 , 3 18 



5 e 

Z. 73 s 10 “ 5 

9 . 8 

US. 1 

0.9534 

0 . 994 

0 , 49 

3.71 

0.0300 

0.087 

0 . 080 



1 j 

a 

P 

4 . 32 X 10' 6 

159.9 

70 . 2 

0.9998 

0.979 

0,88 

3.77 

0. 0082 

0 . 206 

1.043 



^p, c 

0.00290 

5,6 

147.0 

0.9420 

0.937 

0 . 80 ! 

2.87 

, 0. 0717 

0,086 

0 . 314 



Kc 

2.93 x l<f 5 

3.3 

148 . 3 

0 . 9515 

0.078 

0,45 

2 .S 8 

0.0395 

0.085 

0.080 



a' 

©, c 

0 . 001 S 9 

4 . Z 

147 . S 

0 . 9068 

0.994 

0,70 

3,73 

0.0835 

0,086 

0 . 319 



r ■ ' 
°Fe.c 

1 . 54 X 10 " S 

1.0 1 

149 . 0 

0.9197 

0.994 

0.34 

3.69 

0.0409 

0,084 

0 . 081 

15 

zoo* 

° P 

0.000146 

18.4 

90 . 2 

0.9984 

0 . 688 

0.33 

2.23 

0.0075 

0 . 315 

2,008 



° P 

1 . 87 x 10" 6 

7,4 

95.8 

0. 9827 

0. 978 

0.59 

5.53 

0.0388 

0.295 

1 , 83 ! 



0 

0 P 

2 . 20 x 10‘ 8 

25.3 

86. 5 

0.9985 

0 . 971 

0.45 

5.54 

0. 0098 

0 . 330 

3.046 

15 

300 > 

°P 

0.00331 
_ 7 

51.5 

124,0 

0.9990 

0.721 

0.49 

2, 11 

0.0072 

a. 232 

0.910 



°Fe 

1 . 43 x 10 

18,9 1 

HO. 5 

0.9981 

0.634 

0,29 

1, 88 

0.0055 

a. 201 

0 , 203 



V 

4 . 92 x IQ * 5 

13.0 

143.2 

0.9737 

0,989 

0,88 

5,58 

0,0557 

0 . 199 

0,790 



°r. 

6 . Six Id ' 9 

4.3 

147.8 

0.9595 

0 , 989 

0,31 

5.53 

0. 0227 

0 . 191 

0.202 



c # 
P 

1 . 21 x 10 " 7 

75.1 

122.0 

0. 9994 

0 . 982 

0,73 

5,58 

0.0104 

0 . 262 

1 . 616 


Vor the** Jv and E, the calculated and are iero because of the cutoff factor in Eq. (2, 1 Z). 
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